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Abstract 

The fundamental diversity-multiplexing tradeoff (DMT) of the quasi-static fading MIMO Z interference channel 
(ZIC) with channel state information at the transmitters (CSIT) and arbitrary number of antennas at each node is 
derived. A short-term average power constraint is assumed. It is shown that a simple Han-Kobayashi |T| coding scheme, 
where the 2nd transmitter's signal depends on the channel matrix to the first receiver and the 1st user's transmit signal 
is independent of CSIT, can achieve the full CSIT DMT of the ZIC. We also characterize the achievable DMT of a 
transmission scheme, which does not utilize any CSIT and show that for some range of multiplexing gains, the full 
CSIT DMT of the ZIC can be achieved by it. The size of this range of multiplexing gains depends on the system 
parameters such as the number of antennas at the four nodes (referred to hereafter as "antenna configuration"), signal- 
to-noise ratios (SNR) and interference-to-noise ratio (INR) of the direct links and cross link, respectively. Interestingly, 
for certain special cases such as when the interfered receiver has a larger number of antennas or when the INR is 
stronger than the SNRs, the No-CSIT scheme can achieve the F-CSIT DMT for all multiplexing gains. Thus, under 
these circumstances, the optimal DMT of the MIMO ZIC with F-CSIT is same as the DMT of the corresponding 
ZIC with No-CSIT. For other channel configurations, the DMT achievable by the No-CSIT scheme serves as a lower 
bound to the fundamental No-CSIT DMT of the MIMO ZIC. 

Index Terms 

DMT of MIMO Z interference channel. No CSIT, only CSIR, One sided MIMO IC. 
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I. Introduction 

Significant progress in the direction of a complete understanding of the multiuser interference network has been 
made in the last few years, through the information theoretic performance analysis of various basic building blocks 
of this general network. Different performance metrics, such as the capacity, capacity within a constant number 
of bits. Generalized degrees of freedom (GDoF) and degrees of freedom (DoF), has been used to evaluate the 
performance of these basic channel models. While in terms of accuracy capacity is the best performance metric, 
it is also the hardest one to analyze. For instance, the capacity region of the most basic building block, the 2-user 
interference channel is known only for some very special cases IS, [31, ||4l||5l, ||6l, Q, lH, depending on the 
values taken by the various channel parameters. 

The model used for the channel coefficients of a communication system is another important design parameter 
Most of the capacity characterizations mentioned above or approximate capacity characterizations ||9l, ITOl assume 
time-invariant channel coefficients, whereas the channel coefficients of a practical wireless network vary with time 
and should be modeled as a fading one. Most of the works which assume a fading channel model use a much 
coarser performance metric such as the DoF or the GDoF. These metrics can only characterize the rate scaling 
factor/s with average SNR, of the corresponding channel and do not reveal any information about the reliability 
of communication. Diversity-multiplexing tradeoff introduced in ifTTI for the point-to-point channel captures this 
relation between the rate and the reliability of the communication. Encouraged by the importance of knowing the 
best achievable reliability on a channel while communicating at a particular rate, in this paper we choose DMT 
as our performance metric. Further, as a first step towards understanding the general interference network, we 
choose the 2-user Z interference channel as our channel model in this paper In the 2-user Z interference channel 
(ZIC), 2 transmitters communicate to their corresponding receivers via the same signal space, while only one of 
the transmitters interfere with the other receiver. 

Besides being one of the basic building blocks of the general interference network, ZICs also emerge as the natural 
information theoretic model for various practical wireless communication scenario such as femto-cells lfT2l . Also, 
the ZIC is a special case of the 2-user IC. Thus optimal (with respect to some metric) coding and decoding schemes 
on a ZIC might reveal useful insights for the 2-user IC also. For instance, the optimal DMT of the ZIC is an upper 
bound for the 2-user MIMO IC (with and without CSIT). These facts make the analysis of the ZIC an important 
step towards a better understanding of the general multiuser wireless system. Motivated by the aforementioned facts 
in this paper we analyze the DMT of the MIMO ZIC. However, unlike the DMT framework in a point-to-point 
(FTP) channel, ifTTl where there is a single communication link which can be characterized by a single SNR, in 
a multiuser setting such as the one at hand, it is only natural to allow the SNRs and INRs of different links to 
vary with different exponentials with respect to a nominal SNR, denoted as p. This technique was first used in ||9l 
to analyze the DoF region which the authors referred to as the Generalized DoF (GDoF) region, of the 2-user 
SISO IC. Later, this technique was extended to the DMT scenario of the SISO IC in |[l3l and IH. Following 
similar approach, we allow the different INR and SNRs at the receivers to vary exponentially with respect to p 
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with different scaling factors. We refer to the corresponding DMT as the generalized DTM (GDMT) to distinguish 
it from the case when SNR=INR in all the links. 

In this paper, we first derive the DMT of the MIMO ZIC with CSIT and arbitrary number of antennas at each 
node. The achievability is based on a simple Han-Kobayashi IH coding scheme, where the signal to be transmitted 
by the 2nd user depends only on the channel matrix to the first receiver, whereas the transmitted signal of the first 
user does not use any CSIT. The converse is proved by deriving a set of upper bounds to the achievable DMT, 
from a set of upper bounds to the capacity region of the ZIC. The set of upper bounds to the capacity region 
in turn is obtained assuming a genie aided interfered receiver. The computation of the DMT of the MIMO ZIC 
involves the asymptotic joint eigenvalue distribution of two specially correlated random Wishart matrices which 
was recently derived by the authors in jl5) in a different context. Using this distribution result, the fundamental 
DMT of the MIMO ZIC channel with CSIT is established as the solution of a convex optimization problem. While 
it is argued that in general the optimization problem can be solved using numerical methods, closed-form solutions 
are computed for several special cases. Secondly, we characterize the achievable DMT of a transmission scheme 
which does not utilize CSIT. Comparing the achievable DMT of this scheme with the F-CSIT DMT of the ZIC, we 
identify two classes of MIMO ZICs on which the No-CSIT scheme can achieve the F-CSIT DMT. The first class of 
ZICs have equal number of antennas at all the nodes and a stronger INR than a certain threshold (e.g.,Theorem|6| 
and the second class of ZICs have a larger number of antennas at the interfered node than a certain threshold 
(e.g., Theorem|7]i. The above result thus effectively characterizes the DMT of these channels without CSIT because 
the F-CSIT DMT of the channel is an upper bound to the No-CSIT DMT. For other channel configurations, this 
achievable DMT represents a lower bound to the fundamental No-CSIT DMT of the channel. 

An early work in this direction is lfT6l . where the authors derive an achievable DMT on a SISO ZIC with No- 
CSIT. The DMT of the SISO ZIC with F-CSIT can be obtained from 113, where the DMT (F-CSIT) of the 2-user 
SISO IC was derived. In this work, we focus on the MIMO case. In ifTSl . an upper bound to the DMT of a 2-user 
MIMO IC with CSIT was derived for the case in which all nodes have same number of antennas and the direct 
and cross links have the same SNRs and INRs, respectively. This result, if specialized for the MIMO ZIC, provides 
only an upper bound. Our result will prove that for the special case considered in |18|, this upper bound is actually 
tight on a ZIC with F-CSIT. However, the result of this paper on MIMO ZIC is much more general, in the sense 
that we consider arbitrary number of antennas at each node and arbitrary scaling parameters for the different SNRs 
and the INR of the system. Moreover, we also characterize the DMT of the channel with no CSI at the transmitters 
for some specific channel configurations. 

Notations: We denote the conjugate transpose of the matrix A as and its determinant as \A\. C and M 
represent the field of complex and real numbers, respectively. The set of real numbers {x £ M. : a < x < b} will 
be denoted by [a, b]. Furthermore, {x A y), {x V y) and represent the minimum of x and y, the maximum of 
X and y, and the maximum of x and 0, respectively. All the logarithms in this paper are with base 2. We denote 
the distribution of a complex circularly symmetric Gaussian random vector with zero mean and covariance matrix 
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Q as CAf{0,Q). Any two functions f{p) and g{p) of p, where p is the signal to noise ratio (SNR) defined later, 

og(/(p)) _ i;™ log(g(p! 
log(p) - inilp^oo log(p) 



are said to be exponentially equal and denoted as f{p)==g{p) if, linip^oo ^"o^'^f^^"''' = linip-i-oo ^^^^^^ry-- The 



same is true for > and <. 

II. Channel Model and Preliminaries 

We consider a MIMO ZIC as shown in Figure[T] where user 1 (Txi) and user 2 {Tx2) have Mi and M2 antennas 
and receiver 1 (Rxi) and 2 {RX2) have iVi and antennas, respectively. This channel will be referred hereafter 
as a (Ml, A^i, M2, N2) TIC. A slow fading Rayleigh distributed channel model is considered where Hij G ^^j-KMi 
represents the channel matrix between Txi and Rxj. It is assumed that Hn, H21 and H22 are mutually independent 
and contain mutually independent and identically distributed (i.i.d.) CJ\f{0, 1) entries. These channel matrices remain 
fixed for a particular fade duration of the channel and changes in an i.i.d. fashion, in the next. Perfect channel 
state information is assumed at both the receivers (CSIR) and both the transmitters (CSIT). Following lfT4l . we also 
incorporate a real- valued attenuation factor, denoted as rjij, for the signal transmitted from Txi to the receiver Rxj. 
At time t, Txi chooses a vector Xn G C^^' ^ ^ and sends ^/PiXn over the channel, where for the input signals we 
assume the following short term average power constraint: 



N+N(k~l) 

- J2 ""(Q^*) < 1, V fc > 1, z = 1, 2, where Q,t = E (x.tX/,) , (1) 

t=l + W(fc-l) 

where TV represents the number of channel uses for which the channel matrices remain fixed, or in other words the 
fade duration. 

Remark 1: Note that since the transmitters are not allowed to allocate power across different fades of the channel, 
the channel is still in the outage setting, i.e., the delay limited capacity of each of the links of the ZIC is zero |fT9l . 
Thus, the DMT characterization of this outage limited channel make sense. 



Txi(Mi) >; Rxi(Ni) 




RX2(N2) 



Fig. 1: Channel model for the ZIC. 

The received signals at time t can be written as 

Yit = mi^/PiHiiXit + V2i\/P2H2iX2t + Zit, 

Y2t = V22\/l^H22X2t + Z2t, 
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where Z,t € C^'^^ are i.i.d as CN{Q,In,) across i and t. The above equations can be equivalently written in the 
following form. 



Yu = \/SNRni?n^it + VlNR^^2i^2t + Zu] (2) 



Y2t = v/SNR22i/22^2t + Z2t , (3) 

where the normalized inputs XjS satisfy equation ([T]i with equality and SNR^^ and INR^i are the signal-to-noise 

ratio and interference-to-noise ratio, respectively at receiver i, which from now onwards will be denoted by pn and 

Pji, respectively. The performance of a link depends on the strength of the interference relative to the desired signal 

level and when interference strength is much less than the signal strength a better performance. This variation of 

performance due to relative difference in strengths of SNRs and INRs can not be captured through the DMT metric, 

if they differ by only a constant factor. To characterize the DMT under such a scenario we thus let the different 

SNRs and INRs to vary exponentially with respect to a nominal SNR, p with different scaling factors as follows: 

log(SNRn) log(SNR22) 
ail = — — -J—. — , "22 = — — 1—. — , (4) 
log(p) log(p) 

log(INR2i) 

"21 = f-\ ■ 

For brevity, in the sequel we shall use the following notations: SNR^i = pa, INR21 = P21, H = {Hu, H21, H22}, 
P = [pii, /O21, P22] and a [an, a2i, a22]- 

Diversity order of a point-to-point channel ifTTl is defined as the negative SNR exponent of the average probability 
of error at the receiver Since in the present channel model there are 2 receivers and more than one SNR and 
INR parameters, in the next subsection, we provide the definitions of diversity order and the multiplexing gains 
appropriate to the ZIC. 

A. Definition of DMT of a MIMO ZIC 

Let us start by defining the DMT of the channel formally. Let Txi is transmitting information using a codebook 
Ci{p) having 2^^''^''^ codewords, each of length L, at a rate Ri{p), then the corresponding multiplexing gain is 
denoted by where 

T,= lim -^^4' fo"" * = 1'2- (6) 

p-i-oo lOg(p) 

Remark 2: Note that the maximum asymptotic rate supportable by the first direct link is R^^^ = min{7\/i, A^i} log(pii), 
when the second user is silent. Putting this into equation (|6]l we have 

n = min{Mi, A^ijan > min{Mi, A^i}, if an > 1. (7) 

Apparently, it seems that the direct hnk can support a multiplexing-gain strictly larger than min{Mi, TVi}. However, 
this is only a consequence of the fact that the multiplexing gain ri in ^ is defined with respect to (w.rt.) the 
nominal SNR p and therefore, a by product of the more general mathematical model that we assume in here in 
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this paper. With respect to the direct Hnk's SNR pu the multiplexing-gain is still min{A/i, A'^i}, irrespective of the 
value of ail, i-e., 

f, = lim ,^]^^\ = min{Mi,iVi}. (8) 

Pll-^oo l0g(/3ii) 

Alternatively, this apparent difference can also be removed by equating the nominal SNR p to the direct links SNR 
pii which amounts to setting an = 1 in equation (|7]i. 

Now, to define the diversity order, let 'Pe,Ci{p,ri,T2) represents the maximum of the average probability of 
errors at the receivers (averaged over the random channel, Gaussian additive noise at the receivers and different 
codewords of a codebook) at a multiplexing gain pair (ri,r2) and SNR of p, and V*{p,ri,r2) represents the 
minimum 'Pe.Ci{p,i'iii"2) among all possible coding schemes, i.e., 

7'*(p,ri,r2) = min "Pe.cCp, ^'i, ^'2), (9) 

{All possible coding scheme, Ci (p) } 

then the corresponding diversity order ifTTl is defined as 

/ N -log {r*{p,ri,r2)) 
dz\c{Ti.r2)= lim — . (10) 

p^oo 10g(pj 

Note that the diversity order <^zic(''ii ''2) is a function of the relative scaling parameters of the different links a. 
However, for brevity of notation, we shall not mention them explicitly. 

The typical approach to characterize the DMT of a communication channel, whose exact instantaneous end-to-end 
mutual information (IMI) is not known, is to find an upper and a lower bound to it. Then, from this upper and 
lower bound to the IMI a lower and upper bound to the appropriately defined outage event is derived, respectively. 
If the later set of bounds have identical negative SNR exponents then that represents the DMT of the corresponding 
channel. In this paper we adopt the same approach and therefore, need a subset and a superset to the IMI region 
of the channel, which we specify in the next subsection. 

B. A Subset and a Superset to the instantaneous mutual information region 

In this subsection, we shall first derive a set of upper bounds to the various end-to-end mutual information 
defining a superset to the IMI region of the MIMO ZIC. Next, we shall propose a simple superposition coding 
scheme, which can achieve a IMI region, with its various bounds within constant (independent of SNR and channel 
coefficients) number of bits to those of the superset. These bounds will then be used to derive the fundamental 
DMT of the channel. 

Lemma 1: The IMI region of the 2-user MIMO ZIC with F-CSIT, for a given realization of channel matrices 
H, denoted by C('H, p), is contained in the set of real-tuples TZ^{H, p), where 72."('H, p) represents the set of rate 
pairs {Ri,R2) such that Ri,R2>0 and satisfy the following constraints: 

'In^+PuHuH}^ ^e {1,2}; 



i?, < log 
i?i + i?2 < loe 



In^ + p2iH2iHl^ + piiH II 



log 



+ P22H22 {Im2 + P21HI1H21 ) 



bs- 
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In what follows, we find a subset to the IMI region of the channel. Consider a coding scheme where the first 
transmitter uses a random Gaussian code book and the second user uses a superposition code as follows 



(11) 



where U2 (hereafter mentioned as the private part of the message) and W2 (public part of the message) are mutually 
independent complex Gaussian random vectors with covariance matrices as follows: 



EiX^xl) = Im,, E{W2W^) - ^ and 



IE(C/2C/t) 



'AI2 



H21H2I 



(12) 



Remark 3: Note that this covariance split satisfies the power constraint in equation ([T]l. 

Remark 4: The above described coding scheme is clearly a special case of the Han-Kobayashi coding scheme 
where the first transmitter's message does not have any private part. Also the DMT characterized in this paper 
represents the best DMT achievable on a ZIC when both transmitters have full CSIT. However, it will be shown 
shortly that the above coding scheme which use the knowledge of only H21 at Tx2 can achieve the F-CSIT DMT. 

Lemma 2: For a given channel realization Ti, the above described coding scheme can achieve an IMI region 
72.' ("H, p), where 72.' {H, p) represents the set of rate pairs R2) such that i?2 > and satisfies the following 
constraints: 



Ni 



R, < log 
Ri +R2< log 
log 



2iV,; 



lu, «e {1,2}; 

t 



-1 



PuHuHI 
Ini 

+ P22H22 i^hh + P2iH\iH2ij H. 
~2{m+N2)^Iis. 

It has been explained earlier that the superset, TZ^{H, p) given by Lemma [T] and the subset, 72' (7^, p), given by 
Lemma |2] to C{7i,p) can be used to derive a set of upper and lower bounds to DMT of the channel, respectively. 
In what follows, we prove this fact formally. Let Txi is operating at a rate Ri{p) bits per channel use where the 
corresponding multiplexing-gain is = limp^oo for * = 1, 2 and the outage event O is defined as follows 



0^{H:iRi,R2)(^C{H,p)}. 
Following similar method as in ifTTI it can be easily proved that 

V:{p,n,r2)^PT{0), 



(13) 



(14) 
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where 7^* {p,ri,r2) represents the minimum average probabiHty of error achievable on the ZIC, as defined in 
subsection III-AI Now, from Lemma [T] and [2] for any realization of the channel matrices H we have, 

7^'(H,p)CC(H,p)C7^"(H,p); 
Or, {n : (i?i,i?2) i C O C : (i?i,i?2) ^ n\n,p)}; 

Or, Pr{(i?i,i?2) ^7^"(H,p)}<7':(p,rl,r2)<Pr{(i?l,i^2) ^7^'(H,p)}; 
Or, Pr i?2) ^ 7^"(H, p)} <p-'^-c(-i-2)< pr |(^^^ ^2) ^ 7^'(H, p)} ; 
Or, Pr{U,{4, < i?.}}<p"''^'^''"^^'-^)<Pr{U,{/H < R^}} ; 
Or, max Pr {4, < i?J <p-''^'c(-i.'-2)< y Pt{Iu<R,}, 



(15) 



j=l,2.s 



Or, max Pr {/b^ < i?.,} <p"'^z*ic('-i,'-2)< j^ax Pr{/;, <i?,}, (16) 
te{i.2,s} ie{i,2.s} 

where Rg = + R2) and /bj's and /ji's are as defined in Lemma [T] and Lemma |2] respectively. Note that, 

hi = hi + Ci, where Ci is a constant independent of p, for all i E {l,2,s}, which becomes insignificant at 

asymptotic SNR. Therefore, at asymptotic values of p equation ( fTSI l is equivalent to 



max Pr{/b,, < Ri} , 

ie{l,2,s} 



which can be written as 



diiciri,r2) =mindoiiri), where (17) 

i^X 

aOi(?-j )= lim ; — , (18) 

for alH G I = {1, 2, s} and = (ri + r2). 

The only remaining step to characterize the DMT completely is to evaluate the probabilities in equation ( fTSI ). 
which in turn requires the statistics of the mutual information terms /;,/s. It will be shown in the next section 
that this statistics and therefrom the DMT of the channel can be characterized, if only the joint distribution of the 
eigenvalues of 2 mutually correlated random Wishart matrices are known. 

III. Explicit DMT OF THE ZIC 

In this section, we shall evaluate the different SNR exponents, (io;(ri)'s, of the various outage events given in 
equation dTTI l. which would yield the explicit DMT expression for the ZIC. Substituting the right hand sides of the 
first and second bound's in Lemma [T] in equations ( fTSl l we have 



doAn)^ lim ^ I e {1,2}, (19) 

where Ui^fc's are the negative SNR exponents of the ordered eigenvalues of the matrices HuHJ^. The joint distribution 
of {yi,k}'k=i^^"^'^ was specified in ifTTI . Using this distribution and a similar technique as in ifTTI . it can be shown 
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that 

min{A-fi,A''i} 

do, in) = min ^ {M + N + I - 2k)v,^k (20a) 

k=l 
min{A/i,JVi} 

subject to: ^ {a^i - ■Lii,fc)+ < r; (20b) 

< Vi^i < ■ ■ ■ < Vi,-^in{Mi.Ni}- (20c) 

Similar optimization problem will recur in the remaining part of the paper. So we formally state the solution of the 
above problem in the following Lemma for convenience of reference later. 

Lemma 3: If d{r) represents the optimal solution of the optimization problem, 



then. 



mm^{M + N +1 -2i)fi, (21a) 
1=1 

m 

subject to: - < r; (21b) 

1=1 

< Ml < ■ ■ • < Aim, (21c) 
d{r) ~ adM,N ^— ) , for < r < ma, (22) 



where m = min{A/, and dM.N{i") represents the DMT of a. M y. N point-to-point channel and is a piecewise 
linear curve joining the points [M — k){N — k) for k ~ 0,1, ■ ■ ■ m. 
Proof: Putting /i^ = ^ in the optimization problem ( |2TI ) we get, 

d(r) 



: mm 



i=i 



m^{M + N +l~2i)fi[ (23a) 



subject to: "Vfl - ^'M < -; (23b) 
^-^ a 

1=1 

0<m'i < ••• <M™- (23c) 

The solution of this modified optimization problem was derived in [Tl] and is given by 

d(r) , / r\ r 
— =dM.N - , for < - < m, 
a \ a/ a 



or, d{r) =adM,N (^—^ , for < r < ma. 



The solution of the optimization problem ( l20b is now evident from Lemma [3] and is given by 



do. in) = audAu,N. — , V r, e [0, min{M„ N,}au] and i e {1, 2}, (24) 
\aii J 

where dm,nir) is the optimal diversity order of a point-to-point MIMO channel with m transmit and n receive 
antennas, at integer values of r and is point wise linear between integer values of r. To evaluate do^ (fs), we write 
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the bound Its of Lemma [T] in the following way 

hs = log 
+ log 



(^In2 + P22H22 [hh + P2iHI^H2i^ ^ Hi 

log (/tVi +P2li^21^^1i 



Hu 



(a) 



12 



^(1 + A.) + ^(1 + + E(i + p"'"'^^) 



i=l .7 = 1 fc=l 

where in step (a), p = min{M2, A''i}, qi = min{7\/i, iVi}, q2 ~ min{A/2, -^'^2} and we have denoted the ordered 
non-zero (with probability 1) eigenvalues of M^i = {Jni + P2iH2iH\^ Hu, W2 = H22 (^Ihh + P21HI1H21 
and W3 = i?2i^^2i by Ml ^ ■ ■ ■ 5^ Pgi > 0, tti > • • • > 7rg2 > and Ai > • • • > Ap > 0, respectively. Now, using 
the transformations Ai = p^'"', for I < i < p, pj = p^^^, for I < j < qi and TTk = p 1 < fc < 92 in the 
above equation and substituting that in turn in equation ( fTSl l we get 



91 



^21 



i=i 



92 

E( 

A;=l 



^22 - Ik ) 



(25) 



To evaluate this expression we need to derive the joint distribution of 7, f3 and v where 7 = {71 , • • • , 7g2 } and 
similarly v ~ {vi, • • • , Vp} and /3 = {/3i, • • • , (3q-^}. Note that Wi, W2 and are not independent and hence 
neither are 7, /3 and u. However, this distribution can be computed using Theorems 1 and 2 of ifTSl . Using this joint 
distribution, equation dZST l and a similar argument as in IfTTI do^{rs) can be derived as the solution of an convex 
optimization problem as stated in the following Lemma. 

Lemma 4: The negative SNR exponent of the outage event corresponding to the sum bound in Lemma |2] i.e.. 
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dosi'^'s), is equal to the minimum of the following objective function: 

P Ql 

i=i j=i 

92 

+ E(A/2 + A^2 + 1 - 2fc)7fc - {Ml + N2)pa2i 

k=l 

92 rain{{M2-k),N2} qi min{{Ni-j),AIi} 

+ £ {a2i-v^--fk)+ +J2 («2i - w«-ft)+; (26a) 

A:=l j=l i=l 

P 91 92 

constrained by: E(a2i - v^)"^ + ^(aii - ft )''' + ^("22 - Jk)^ < r^; (26b) 

i=l j=l k=l 

0<vi<---<Vp; (26c) 

0<f3i<---<P,,; (26d) 

0<7i < ••• < 792; (26e) 

(u, + ft-) > a2i, V(i + j) > (iVi + 1); (26f) 

(w,: + 7fc) > "21, V(^ + fc) > (A/2 + 1). (26g) 

Proof: The proof is relegated to Appendix |A] ■ 
Differentiating the objective function with respect to {vj}, {ft } and {7fe}, it can be easily verified that ( I26al ) is 
a convex function of these variables. The constraints on the other hand are all linear. Therefore, equation ( l26b 
represents an convex optimization problem (e.g., see subsection 4.2.1 in ll20l ) and hence can be solved efficiently 
using numerical methods. Since we have already found expressions for doi for i = 1, 2 as in equation (l24l i. Lemma|4] 
provides the last piece of the puzzle required to characterize the DMT of the MIMO ZIC, by evaluating rfo, • This 
is stated formally in the following Theorem. 

Theorem 1: The DMT of the (Afi, A^i, M2, A^2) ZIC, with FCSIT and short term power allocation scheme, is 
given as 

dzic^^^ (n , r2 ) = mill do,{ri), 

'iG|1.2,s} 

where do. (r;) for i = 1,2 and i = s are given by equation ( l24b and Lemma ID respectively. 

Although the computation of da^ (r,, ) and hence characterization of the DMT of a general ZIC with arbitrary 
number of antennas at each node require application of numerical methods, in what follows, we shall provide 
closed form expressions for it for various special cases. Since the DMT with FCSIT acts as an upper bound to the 
DMT of the channel with no-CSIT or only CSIR, these expressions facilitates an easy characterization of the gap 
between perfect and no CSIT performances of the channel. The no-CSIT DMT of the channel will be characterized 
in section UyI for different range of values of 021 and no of antennas at different nodes. 

The central idea in all of the following proof is the fact that the steepest descent method provides a global optimal 
value of a convex optimization problem: it is well known that the steepest descent method provides a local optimal 
value of the objective function. However, a local optimal solution is equal to a global one as well, for a convex 
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function 1201 . Combining the above two facts we conclude that the value obtained by the steepest descent method 
is actually the global minimum of the objective function. The first case considered is a class of channels where all 
the nodes have equal number of antennas. 

Theorem 2: Consider the MIMO ZIC with Mj ^ M2 ^ Nx = f^2 = n and SNRs and INRs of different hnks are 
as described in Section HIl with a^i ~ ct and 022 = 1 = an- The best achievable diversity order on this channel, 
with F-CSIT and short term average power constraint ([T]i, at multiplexing gain pair (ri,r2), is given by 

4icTn^n,„,«)(''i'''2) = min{rfn,„(?'i),rf„,«(?^2),<^(?J^„,„,„)(r-.)} (27) 

where do^rg) for the special channel configuration being considered is denoted by d^J^^^^^ ^ if " < 1, 

then 

,FCSiT r N J (^dn,3n{^) + 2«'(1 " «), for < < na; 



I 2(1 - a)d„,n( 2(i-a) )^ for na < < n{2 - a). 



and if 1 < a, then 



dn -^n (r,)+'n'^(a^l), < r, < n: 
jFCSIT I'^^-j ^ J "^•^"^ ' ^ J, - s - , ^2^^ 

\ (a - l)d„_„ (2^) , n < < na. 
Proof of Theorem |2} The proof is relegated to Appendix |B] ■ 

Remark 5: Note that for a = 1, the optimal DMT becomes d^^(^'^_„_^j(rs) = min {d„,„(ri), (i„,„(r2), d„,3„(ri +r2)} 
which is exactly the upper bound derived in ifTSl . 



A. The DMT of a Femto-Cell 

A practical communication channel following the ZIC signal model appears in the so called Femtocell envi- 
ronment. The Femtocell |T2| concept is an outcome of the telecommunication industry's efforts to provide high- 
throughput, high quality services into the user's home. Consider the scenario depicted in figure |2] where the larger 
circle represents the macro cell serviced by the macro cell base station (MCBS). Within this macro cell is the 
smaller circle represents a small area where the signal from the MCBS either does not reach with enough strength 
or does not reach at all, hereafter referred to as the Femtocell. To provide coverage in this region a smaller user 
deployed base station connected to the backbone can be used, which is called the Femto cell BS (FCBS). This FCBS 
can provide mobile services to the users within the Femtocell just like a WiFi access point. The basic difference 
between the FCBS and WiFi access point is that the former operates in a licensed band. 

Now, let us consider the downlink communication on such a channel with one mobile user in both the Femtocell 
and the macro-cell. Note that since the MCBS signal does not reach the mobile user within the Femtocell, the 
signal input-output follows the ZIC model. To model the larger SNR of the Femtocell direct link we can assume 
that 0:22 > 1. In what follows, we shall derive the DMT of this channel. 

Theorem 3: Consider a ZIC with Mi = M2 = Ni = N2 = n, 022 = a > 1 and an = 021 = 1, F-CSIT and 
short term average power constraint given by The optimal diversity order of this channel at a multiplexing gain 
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pair (ri,r2) is given by 



where do^{rs) for the special channel configuration being considered, is denoted by <is™'°„ „ „)(''«), and 

,Femto ( X I dnMrs) + ^^(q! - 1), for < < n; 

y [a - l)a«,K( ), for n < < na. 

Proof: The proof is relegated to Appendix |C] ■ 
Remark 6: Note that the fundamental DMT of the ZIC with single antenna nodes and a22 = ct, ctii = ct2i ~ 1 
was derived in |2T|. This clearly is a special case of Theorem [3] and can be obtained by putting n = 1. 

Typically, in a multiuser communication scenario one end - say the base station in a cellular network - can host 
more antennas than the other, in what follows, we consider a case where Mi = M2 = M < min{A^i, A^2}- 

Theorem 4: Consider the ZIC with AIi = 1\I2 = M < minjA'^i, N2}, an = 0122 = ci2i = 1, F-CSIT and short 
term average power constraint given by ([T]i. The optimal diversity order achievable on this channel at a multiplexing 
gain pair (ri,r2) is given by 

C^M^/Vi,M,JV2(''l:'"2) = min|dM,Ari(ri),dM,JV2(?'2),d^^(^7Vi3/,W2) 

where doX^s) for the special channel configuration being considered, is denoted by d^J^^^J 1/^2)^^'*)' ^^'^ 

dAUM+N,+N,)irs) + M{Ni - M); < < M; 
d2M,N^irs); M<r,< inm{Ni,2M}. 
Proof: The proof is relegated to Appendix iP] ■ 
Next we shall demonstrate the performance loss in terms of DMT, due to lack of CSIT. In Figure [3] explicit F- 
CSIT DMT curves for a few antenna configurations are plotted and compared against the performance of orthogonal 
schemes such as frequency division (FD) or time division (TD) multiple-access which do not require CSIT. It can 



jFCSIT / x 



December 6, 201 1 



DRAFT 



14 




Symmetric Multiplexing gain, r_ 



Fig. 3: Optimal DMT of different ZICs with a = [1, 1, 1]. 



be noticed from the figure that the gain due to CSIT, over the orthogonal access schemes can be significant, 
particularly in MIMO ZICs. While using better coding-decoding schemes this gap can be reduced, in general, with 
CSIT a better performance can be achieved. However, to evaluate this gap in performance due to lack of CSIT, 
exactly it is necessary to know the best DMT achievable on the channel without any CSIT. Popular approaches 
of characterization of the No-CSIT DMT involves either evaluating the optimal instantaneous mutual information 
region of the channel without CSIT exactly or finding subsets and supersets of it which are within a constant 
number of bits to each other, as was done in the case of F-CSIT considered so far. For the no CSIT case either of 
the above two objectives are very hard to achieve. However, in the following section, we bypass this approach and 
charcterize the No-CSIT DMT of the ZIC for some specific values of the different channel parameters such as the 
number of antennas at the different nodes and the INR parameter a. 

IV. DMT WITH No CSIT 

To avoid the above mentioned difficulties in this section, we use as the F-CSIT DMT derived in the previous 
sections, an upper bound to the No-CSIT DMT of the channel. Then, we derive the achievable DMT of a No-CSIT 
transmit-receive scheme, which for two special classes of ZICs meets the upper bound and therefore, represents the 
fundamental No-CSIT DMT of the corresponding ZICs. In what follows, we shall describe the No-CSIT transmit- 
receive scheme first. 

Let both the users encode their messages using independent Gaussian signals. On the other hand, since Rxi is not 
interested in the signal transmitted by Tx2 we consider a decoder at Rxi which does joint maximum-likelihood (ML) 
decoding of both the messages. However, the event where only the second user's message is decoded incorrectly 
is not considered as an error event. Rx2 uses an ML decoder to decode its own message. Hereafter, we will refer 
to this scheme as the Individual ML (IML) decoder and the encoding-decoding scheme as the Independent coding 
IML decoding scheme or the IIML scheme. 
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An achievable rate region of the individual ML decoder and independent Gaussian coding at both the users is 
given by the following set of rate tuples 



"^IML — 



{(i?i,i?2) -Ri < logdct (in, + JfH^iHl^ = Ic^; (30a) 



Ml 

_P_ 

M2 



i?2 < log det ( In, + -^H22HI^ ) ^ ; (30b) 



+ ra) < logdct [In, + Jf^HnHi^ + ^^H2iHlj ^ 1^/,] (30c) 



Using the above expression for the achievable rate region, the corresponding achievable DMT of this transmit- 
receive scheme, i.e., mutually independent Gaussian coding at each transmitters and IML decoder at the interfered 
receiver, can be easily computed using standard techniques. The result is specified in the following Lemma. 

Lemma 5: If we denote the achievable diversity order of the IML scheme, at multiplexing gain pair (ri,r2), by 
^(^';,iVi,M,,iv.)(^i'^2) then. 



d'(t\,,N,,M.,N,M^r2) > . mill ^,^_^^)(r,)}, (31) 

where Tg = (ri + r2) and 

jIML / \ r log (Pr (/ci ^ fi)) w ■ ^ n o 1 /TON 

d^X^h,N,.M■2,N,)ir^) = ^lim ^—^^ , V i e {l,2,s}. (32) 

Proof: The proof is relegated to Appendix |E] ■ 
Note that 1^ and /^^ represents the mutual information of a point-to-point channel with channel matrices Hn and 
H22, respectively. Therefore, by the results of ifTTI we have 

<f^^h,N^,M2.N,)ir^) = du^.N^). < r, < min{M„iV,}, (33) 

and i e {1, 2}. To analyze the outage event due to the third bound of the achievable rate region, we first approximate 
/c^ by another term which does not differ from it by more than a constant. Note that 



logdct [In, + pHuHl + p"H2iHl,j - N, log(max{Mi, Afa}), 

< Ic^ < log dct (/at, + pHuHl^ + p"H2iHI) = I',^ . (34) 

Since a constant independent of p, does not matter in the high SNR analysis, to compute we can use in 
place of ■ Next, we write in the following manner: 

-log dct (/at, + pHnHl, + p"i/2ii^li) , (35) 

-log det (/Ml + pHl^Hii) + log dct (In, + p"ff2ii/]i) , (36) 



(a) 



5](l+p""Ai,) + ^(l + p"-A,)[, (37) 

]=i i=i 



where Hn = [Ini + p"-f^2i^^2i j " ™d in step (a), p = miii{Ajf2, -/Vi}, qi = min{i\/i A^i}. Also, we 
have denoted the ordered non-zero (with probability 1) eigenvalues of Wi — \ Ini + P21H21III1) Hii and 
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W2 = H21H21 by /.Ji > • • • > iiq^ > and Ai > • • • > Ap > 0, respectively. Now, using the transformations 
Xi = p^^'\ for 1 < i < p, ^ p^^\ for 1 < J < 9i in the above equation and substituting that in turn in equation 
( |32] | we get 

p-d™(Mi,«i,M„«,,('-.)^Pr 1^1^(021 +^(an 1 < r.j . (38) 

To evaluate this expression we need to derive the joint distribution of /? and v. Note that, since Wi and W2 are 
mutually correlated and so are (5 and v. As already stated earlier, in general characterizing the joint distribution of 
the eigenvalues of such mutually correlated random matrices is a hard problem. However, in what follows, we show 
that this distribution can be computed using Theorems 1 and 2 of ITSl . which in turn facilitates the characterization 
of ^(r,). 

Lemma 6: The negative SNR exponent of the outage event corresponding to the sum bound in ( l30l l. i.e., 
'^^Mi Ni M2 No)^^"")' equal to the minimum of the following objective function; 

p qi 
dTtMuNuAh,N.)irs) = min ^(M2 + N, + M, + 1 - 2z)v, + ^(A/i + N, + 1 - 2j)(3, 

(^'■^) ^=l j = l 

qi min{{Ni~j'),Mi} 

-Mipa2i+J2 ("2i-w»-/Sj)+; (39a) 

j=l i=\ 

p 11 

constrained by: ^(a2i - Vi)^ + ^(aii - f3j)'^ < r^; (39b) 

i=l 3=1 

0<vi<---<Vp; (39c) 

0</3i<...</3,^; (39d) 

(u, + I3j) > a2i, V(^ + j) > (Ni + 1). (39e) 

Proof: The proof is relegated to Appendix IF] ■ 
Theorem 5 (A lower bound to the No-CSIT DMT of the ZIC): 1) The optimal diversity order achievable by 
the IIML scheme described above, on a {Mi, Ni, M2, N2) ZIC without CSIT is given as 

rfLB:zic(n,r2) = . min rfJ^L„Ar„M„Ar,) 

where df^^j^ ^'^^ i = 1, 2 and i = s are given by equation ( [33] ) and Lemma |6l respectively. 

2) This DMT also represents a lower bound to the No-CSIT DMT of the (Mi, TVi, M2, iV2) ZIC. 

Proof: The first part of the Theorem follows from Lemma |5] and the second part of the Theorem follows from 
the fact that the IIML scheme is only one of the numerous transmit-receive schemes that are possible on the ZIC. 

■ 

Although the computation of rf™'^^^ M2 N2)^^^^ ^^'^ hence characterization of the lower bound to the No-CSIT 
DMT of a general ZIC with arbitrary number of antennas at each node require application of numerical methods, 
in what follows, we shall provide closed form expressions for it for various special cases. We shall see that for two 
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special classes of ZICs this lower bound meets the upper bound, i.e., the FCSIT DMT of the channel. We start 
with the case where all the nodes have equal number of antennas. 

Lemma 7: On a MIMO ZIC with n antennas at all the nodes, an = 022 = 1 and a^i = a, the IML decoder 
can achieve the following DMT 



4n!^i,n,n)(^i'^2) = min |d„,„ (ri ) , d„,„ (r2 ) , d™^^„^„_„) (r, ) | 
where d^(^,„,„,„)(r.) is given as 



,iML ( X , dnM^s) + n^ia - 1), < < n; 

ds,{n,n,n,n){T^s) = < (40) 

[a - l)d„,„(^3f)' n<'rs< na. 
Proof: The desired result is obtained by following the same steps as in the 2"'^ part of Theorem |2] ■ 
Figure S] illustrates that the IML decoder can achieve the DMT (with F-CSIT) of the MIMO ZIC on a region of 
low multiplexing gains. 




Comparing equations ( [29b to ( l40b we see that the IML decoder can achieve the F-CSIT DMT for high multiplexing 
gain values also, when a > L This fact raises the natural question: is it possible for the IML decoder to achieve the 
F-CSIT DMT for all multiplexing gains and if it is, under what circumstances? It turns out that if the interference 
is strong enough then the IML decoder can achieve the F-CSIT DMT for all symmetric multiplexing gains. For 
example. Figure |5] illustrates this effect on a ZIC with 2 antennas at all the nodes. 

This characteristics of the DMT on MIMO ZICs for general n is captured by the following Lemma. 

Theorem 6: The DMT specified in Lemma |7] represents the fundamental DMT of the channel with No-CSIT and 
ri ~ r2 ~ r, if 

a>l + ^^t^. (41) 
Proof: Detailed proof will be provided in Appendix |G] ■ 
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0.5 
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0.2 0.4 0.6 0.8 1 1.2 

Symmetric multiplexing gain, ^^=^2'^'' 
Fig. 5: The DMT on the (2, 2, 2, 2) ZIC with different a. 

Lemma 8: Consider the MIMO ZIC as in Theorem |4] but no CSI at the transmitters. The achievable diversity 
order of the IIML scheme on this channel, at a multiplexing gain pair (ri,r2), is given by 

^(A/;Afi,A/,Af2)(''l'''2) = min{dM,JVi(r-l),dj\/,Ar2(''2),d2M,Afi(?'s)}- 

Proof: From Lemma |5] and equation (l33l l it is clear that to prove the Lemma it is sufficient to derive 
an expression for c^^'^.j MAf2)(^'*)' Towards that, for convenience we use I^^ instead of /c, to evaluate the 
corresponding outage event since the two are within a constant number which does not matter at asymptotic SNR 
as was shown in equation ( l34l i. Using this in equation (|32] | along with the facts that Mi = M2 = M and a ~ 1 
we get 

p-<U.N,.M.N,>ir,)^p,. jiog^gt (/at, + pHnHl, + pH2iHl^) < log(p)} , 

=Pr {log det {In, + pH,Hl) < r, log(p) } , (42) 

where He = [Hn H21] G £,Nix{2M) identically distributed as the other channel matrices, since Hn and H21 
are mutually independent. However the right hand side of the last equation represents the outage probability of 
an A^i X 2M point-to-point MIMO channel whose negative SNR exponent was computed in ifTTI and is given by 
rfjVi,2j\f(?'s)- Using this in equation (|42] | we get, 

p~'^™('j\f .JVi .M,N2) ) z^p~'^"l ('"s ) 

or, 4'^(A/,Ari,A/,W2)(''s) ^dN,,2M{rs)- 

Substituting this and equation ( |33] | into equation (IJTT i we obtain the desired result. ■ 
Figure |6] depicts the comparison of the achievable DMT of the IML decoder with that of the fundamental F-CSIT 
DMT of the channel on two different MIMO ZICs. Comparing the performance improvement of the IML decoder 
on the (3, 4, 3, 3) ZIC with respect to that on the (3, 3, 3, 3) ZIC, we realize that a larger number of antennas at 
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Symmetric multiplexing gain, i'^=''2=r 



Fig. 6: Effect of a larger number of antennas at the interfered node on ZICs with a = [1, 1, 1]. 



the interfered receiver can completely compensate for the lack of CSIT. Again by the argument that the FCSIT 
DMT represents an upper bound to the No-CSIT DMT of the channel, the observation from Fig. |6] imply that the 
DMT of Lemma [8] represents the fundamental DMT of the (3, 4, 3, 3) ZIC with only CSIR. It turns out that, this 
channel is only a member of a large class of ZICs for which the No-CSIT DMT can be characterized. This class 
of channels is specified in the next Theorem. 

Theorem 7: The DMT specified in Lemma |8] represents the fundamental DMT of the channel with No-CSIT and 
ri = 7'2 = ^, if 

^ - M 

Proof: Comparing Theorem |4] and Lemma [8] the desired result can be obtained following the similar steps as 
in the proof of Theorem |6l ■ 

V. Conclusion 

The DMT of the MIMO ZIC with CSIT is characterized. It is shown that the knowledge of H21 at the 2nd 
transmitter only is sufficient to achieve the F-CSIT DMT of the channel. The No-CSIT DMT of two special class 
of ZICs have been characterized revealing useful insights about the system such as: a stronger interference or a 
larger number of antennas at the interfered receiver can completely compensate for the lack of CSIT on a ZIC. 
Characterizing the DMT of the general ZIC with no CSIT is an interesting problem for future research. 

Appendix A 
Proof of Lemma|4] 

Recall that, the negative SNR exponents of W^i, W2 and Wz are denoted by {Pj}'^^^, {'^k\'k=i {''^jlfLi' 
respectively. Interestingly, these matrices have exactly the same structure specified in Theorem 1 and 2 of ifTSll . 
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except from the yo"^i which can be easily taken care of. Therefore, using these Theorems we obtain the following 
conditional distributions: 



fwi\WsW\^)^\ - (43) 

if il3,v)iBi 



where 



{91 / min{{Ni~j),Ah} \ ^ ] 

J2 I (^^1 + ^1 + 1 - 2j)/?j + ("21 -V,- /3,)+ 1 - Ml ^(^21 - V,)+ \ ,(44) 

6i={(/3,D): 0<vi<--- <Vp; < A <•••</?,, ; {v, + /3j) > a2i, V(2 + j) > (A^i + 1)} (45) 

and 



fw-i\wAl\^)= { (46) 
if {i,v)iB2 



where 



mm{{M2-k),N2} 

E2{l: v)={Y, \ {M2 + N2 + I- 2k)jk + ("21 - - 7fc)+ I - ^2 5I(«2i - V^)+ \ ,(47) 

k=l \ 1=1 / i=l 



B2 



= |(/3,u) : 0< ui < ••• < Up; < 71 < ••■ <792; (w»+7fc) > "21, V(i + fc) > (M2 + 1)} (48) 



These pdfs are easily obtained from equation (10) of [151] by changing it properly due to the presence of the 
instead of p in flJI and using the fact that 

3-ELi(M2+iVi+i-2zK if 0<ui < ••• <z;„; 
fwM={ - - - P, ^49^ 

otherwise. 

It was proved in 1221 that given the non-zero eigenvalues of W3, Wi and W2 are conditionally independent 

random matrices. Since the non-zero eigenvalues of ^21^^21 '^^'^ h\^H2i are exactly the same for each realization 

of the matrix i?2i, given v, /3 and 7 are conditionally independent of each other. Intuitively, it is a well known fact 

in the literature of random matrix theory that the eigenvalues of Wi are dependent on the matrix W3 only through 

it's eigenvalues. The same is true for W2- Therefore, given the eigenvalues of W3, the eigenvalues of Wi and W2 

are conditionally independent. Using this fact we have the following: 



/wi,W2,W3(/3,7,'y) =fwuW2\W3{l3,Mv)fwA'L'), (50) 
= fwi \W:i(i3\v)fw2\W3i^\v)fw3{.v). (51) 

'For a detailed proof of this fact the reader is refen'ed to Lemma 1 of 1221 . 
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Now, substituting equations ( |43] l. ( |46l ) and ( |49l ) into the above equation we obtain the joint distribution fwi,W2,W3 {P, 7, 
When this joint distribution and the Laplace's method is used following the similar approach as in ifTTl to evaluate 
equation dZSl ) we obtain the following: 

do, (rs) = mill + A^i + 1 - 2i)v, + ^^(A/i + Ni + 1 - 2j)pj + ^{M2 + iVj + 1 - 2k)-fu 

1=1 j = l k=l 

p q2 min{{Al2—k),N2} 

- {Mi+N2)J2i^2i~Vi)+ + J2 £ (a2i-u»-7fc)^ 

i=l fc=l i=l 

gi min{(7Vi-j),A/i} 

+ £ E («2i-t^,-/?j)+; (52a) 
i=i 1=1 

P 91 92 

consti-ained by: ^(a2i - Uj)"'" + ^("ii - l^j)^ + ^(022 - 1k)~^ < Ts] (52b) 

i=l j=l fc=l 

< wi < • • • < Up; (52c) 

0</3i<---</3,,; (52d) 

0<7i < ■•• <792; (52e) 

(u, + Pj) > a2i, V(z + j) > (iVi + 1); (52f) 

{v, + 7fc) > a2i, V(i + fc) > (Af2 + 1). (52g) 

Finally, the desired result follows from the fact that by restricting ai < a2i for all i < p does not change the 
optimal solution of the above problem. This can be proved as follows: suppose the optimal solution have > a2i 
for some i. Now, since the objective function is monotonically decreasing function of ai for all i, substituting 
ttj = 0^21 does not violate any of the constraints but reduces the objective function. But, that means the earlier 
solution was not really the optimal solution. Therefore, in the optimal solution we must have a; < a2i for all 
i < p. However, with this constraint we have 

(021 - Ui)+ {a2i - Vi), \f i <p. 
Substituting this in equation ( I52al i we get the desired result. 



Appendix B 
Proof of Theorem[2] 

The main steps of the proof can be described as follows. First, we simplify the optimization problem in equation 
by putting specific values of the different parameters, such as Mi, Ni and aij as stated in the statement of the 
Theorem, in equation (l26T l. Then we calculate a local minimum of this simplified optimization problem for each 
value of r, using the steepest descent method, which by the previous argument then represents a global minimum. 
The later part of the problem, i.e., the computation of the local minimum will be carried out in two steps: in step 
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one, we consider the case when a < I and then in the second step we consider the remaining case. Let us start by 
deriving the simplified optimization problem. 

Substituting Mi = M2 = Ni = N2 = n, an = 0^22 = 1 and 021 = a in the optimization problem of equation 
we obtain the following: 



min ^(4n + 1 - 2i)v,+ ^{2n + I - 2j)Pj +^(2n + 1 - 2fc)7fc 

(•j7,7,/3,a) 



k=l 



n (n-k) n (n-j) 

-2n'a + ^ (a - w» - 7fc)+ + 1] 1] (" " " ft )^; 
fc=i i=i j=i 1=1 



constrained by: '^^(a 



^(i-/3,)+ + ^(i-7fc)+<^.; 



fc=l 



i=l 

< Vi < ■ ■ ■ < Vn] 

0</3i<-- •</?„; 
< 71 < • • ■ < 7n; 

(w, +/3,)>a, V(i+j) > + 
(w,: +7fc) > a, V(i + fc) > (n+ 1). 



(53a) 

(53b) 

(53c) 
(53d) 
(53e) 
(53f) 
(53g) 



As stated earlier, in what follows, we solve this optimization problem in two steps; in the first step, we assume 
a < 1. 



A. Step 1: (a<l) 

To apply the steepest descent method we first compute the rate of change of the objective function with respect 
to the various parameters. Differentiating the objective function in equation ( |53at we obtain the following: 

do 

= (4n + 1 - 2i), 1 < j < n; (54) 

Qi— O,"' ,cti_i— 0,ai_^i— I,--- ,Qti — 1,71 — 1 

= (2n - 1) < (4n + 1 - 2i), V 1 < i < n; (55) 



dai 



971 



Qfl — 1,71 — 1 

= (2n- 1) < (4n+ 1 - 2i), V 1 < i < n. (56) 

Ql = l,/3l = l 

Note that it is sufficient, to consider the decay of the function with respect to (w.rt.) /3i and 71 only, because 
of the decreasing slope of the objective function with increasing index of /3 and 7 and equation (I53dl i and ( I53el ). 
Therefore, it is clear from the slopes of the function given above that for (i — l)a < rs < ia, the steepest descent 
of the function is along decreasing values of a;, while /3i =71 = 1. Note that for these values of /3i and 71 the last 



two terms of equation ( I53ab vanishes and equations (I53dl i-(53g 1 becomes redundant and as a result the optimization 
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problem of ( |53] l simplifies to the following: 

n 

iin^(4n + l-2i)wj + 2?i2(l-Q;); (57a) 



mm 

1=1 



constrained by: — Wj)^ < (57b) 

1=1 

< < • • • < u„; (57c) 
The solution of this optimization problem follows from Lemma |3] and is given by, 

d,{rs) = adn,zn (^) + 277.2(1 - a), < < na. (58) 

The above solution also imply that for a > na, the optimal solution have ai — 0, Vi, which when substituted 
in equation ( |53T l we obtain the following problem. 

n n 

min^(2n + 1 - 2j)/3j + ^(2n + 1 - 2k)jk - 2n^ (59a) 

]=i k=i 

constrained by: ^(1 - + ^(1 - jk)^ < (r^ - na); (59b) 
j=i k=i 

a</?i < ••• </3„; (59c) 
a<7i<---<7«; (59d) 

(59e) 

Note that the last two summands in the objective function ( |53] ) is zero because of equations ( |59c| i and ( |59d| ). Now, 
from the symmetry of the optimization problem ( |59l ) w.r.t. /3i and 7^, we can assume without loss of generality that 
the optimal solution will have /3.i = ji, Vi. Substituting this and 6i ~ a in equation ( [59l ) we get the following 
equivalent optimization problem: 

n 

min2 ^(2n + 1 - 2j)5j , (60a) 

constrained by: ^(1 - a - (5^)+ < f j ; (60b) 

j=i V / 

Q<5i<---<5n; (60c) 

(60d) 

The solution of this optimization problem follows again from Lemma |3] and is given by 

ds{Ts) =2(1 - a)d„,„ f , < < n(l - a) 

( V — TtOL \ 

=2(1 - a)d„,„ -f^ ^ , na<r,< n(2 - a). (61) 



.2(1 - a) 

Combining equations ( fSST l and (1611 1 we obtain equation ( |28] | of Theorem |2] and we have completed the first step of 
this proof. In what follows we consider the remaining case, when a > 1. 
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B. Step 2: (a>l) 

Differentiating the objective function in equation (I53al ) we obtain the following; 



dai 
dOs 



ai=a V i>2,/3i = l,7i = l 
ai— Q V 2>2,7i— 1 



^71 



(4?! — 1), for ai > (q — 1); 
(2n-l), for < ai < (a - 1); 

(2n- 1), for 1); 
ai=a V i>2,/3i=i | (2n — 2), for < Q!i < (a — 1); 



(62) 



(63) 



5a 1 



= (4n - 3),for a2 > (a - 1). (64) 

a.i—{a — l)^a.i—a V ?>3,/3i = l,7i = l 

Comparing equations ( |62l ) and ( |63] ), we realize that for < < 1, the steepest descent is along the direction of 
decreasing ai, while /3i = 71 = 1. On the other hand, comparing equations ( |62] i. ( |63] l and ( |64l i it is clear that 
beyond = 1, decreasing a2 has the steepest descent than /?i,7i and even ai. In the same way it can be proved 
that for {i — \)<rs< i, the steepest descent of the function is along decreasing values of ai, while /3i = 71 = 1. 
Note that for these values of /3i and 71 the last two terms of equation ( 15 3 at vanishes and equations ( |53dt -(53g 
becomes redundant and as a result the optimization problem of ( |53] l simplifies to the following: 



min ^(4n + 1 - 2i)v, - 2n^{a - 1); 

i=l 
n 

constrained by: ^^(a — Vi)^ < r^; 
1=1 

(a - 1) < wi < • • • < w„; 



(65a) 

(65b) 
(65c) 



Substituting v'^ = — (a — 1) for all i in the above set of equations we obtain the following equivalent optimization 
problem. 



min^(4n + 1 - 2i)vi + r?{a - 1); 

i=l 
n 

constrained by: ^^(1 — v'^'^ < 

i=l 

0<v[<--- <v'^, 



(66a) 

(66b) 
(66c) 



which in turn by Lemma [3] have the following optimal value 

dsirs) = adn.3n (rs) + n^{a - 1), < Ts < n. (67) 
It is clear from this solution that for r s > n, ai < {a — \) for all i, and the optimization problem (|53] | reduces 
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to the following: 



\Os^ mill ^(4n+ 1 - 2i)u,+ ^(2n + 1 - 2j)/3j + ^(2n+ 1 - 2fc)7fc 

[vn^PA) j-^i fc^i 

n (n-fc) „ ("-j) 

-2n^a + ^ ^ (a - _ -^,)+ + ^ ^ (a - - /?,)+; (68a) 

A.— 1 i^l j^l i^l 

n n n 

constrained by: ^(q - 1 - v.,) + ^(1 - + ^(1 - 7^)+ < (r^ - n); (68b) 

1=1 j=l k=l 

< ui < • • • < w„ < (a - 1); (68c) 
< /3i < • • • < /3„; (68d) 

< 71 < • • ■ < 7n; (68e) 

(u, + I3j) > a, V(z + j) >{n + 1); (68f) 
(v, + 7fc) > a, V(z + fc) > (71 + 1). (68g) 



Note the upper bound in equation ( I68cb . which is different from equation ( I53cl i. Also, since we are seeking for 
the minimum value of the objective function which decreases with every f3i and jj, in the optimal solution these 



parameters must take their minimum value, which from equations ( |68fl l and ( 68g 1 is given by 



I3j\min = 7j Imin = « - V^-j + l , (69) 

which along with the ordering among the u^'s, /?j's and 7j's imply that all the terms in the last two summands of 
the objective function are non-negative. Substituting these minimum values and thereby eliminating /3j's and 7j's 
from the optimization problem we obtain the following equivalent optimization problem: 

n 

min^(2n + 1 - 2i)vi, (70a) 

i=l 
n 

constrained by: ^^(a — 1 — Ui) < (?'s ~ f^); (70b) 
< ui < • • • < u„ < (a - 1), (70c) 
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where the simplification of the objective function involves the following algebraic computations 

n n 

Os =^(4n + 1 - 2i)v^ + 2 ^(2n + 1 - 2j){a 

n (n-fc) ri (n-j) 

+ (a - - (a - Vn-k+i)) + Y (o^- - (a - Vn-j+i))] 

k=l 1=1 j=l 1=1 

n n 

=^(4n + 1 - 2i)v, + 2^(2n + 1 - 2j){a - Vn-j+i) ~ 2n^a 
i=i j=i 

n n 71 

+n{n - l)a - ^ 2{n - i)vi - ^(a - Vn-k+i) - - - Vn-j+i); 

i=l fc=l j = l 

n n 

=^(2n + l)v, + 2 ^(n + 1 - j){a ~ Vn-j+i) - n{n + l)a 

i=l 3 = 1 

n n 

=^(2n + 1 - 2i)v^ + 2 ^(n + 1 - j)a - n{n + l)a 
i=i j=i 

n 

=^(271+ 1 - 2i)vi. 

i=l 

The solution of this optimization problem follows again from Lemma [3] and is given by 

/ f — \ 

ds{rs) =(a - l)dn,n -r , < (r^ - n) < n{a - 1) 

f — n \ 

= {a - l)d„,„ 7 — , n<rs< na. (71) 

Combining equations (|67] | and (TtTI i we obtain equation ( |29] | of Theorem |2] 

Appendix C 
Proof of Theorem[3] 

When we substitute the specific values of the different parameters in equation ( |26] |. it reduces to the following 



d'Tn:n,n.n)irs) = min^(4n + 1 - 2tyH + Y.{2n + 1 - 2j)/3, + ^(2n + 1 - 2k)^k 

1=1 j=l k=l 

n (n-k) n (n-j) 

- + E E (1 - - ^^)^ + E E (1 - ^» - (72a) 

n n n 

constrained by: ^(1 - t;j)+ + ^(1 - I3j)+ + ^(a - 7^)+ < r^; (72b) 

i=l j=l fc=l 

< ui < • • • < u„; (72c) 

< /3i < • • • < /3„; (72d) 

< 71 < • • • < 7n; (72e) 

iv^+|3J)>l, y{i+j)>{n + l); (72f) 

(t^z+7fc) > 1, V(^ + fc) > (n + 1). (72g) 
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Differentiating the objective function with respect to Vi, V and 71 we find that 

^'^Femto 



dv,. 



= (4n + 1 — 2«) > 

-1,71— a 



= (2n-i), (73) 

71— a,?jfc— 0, yk<.i,Vk—'i-yk>i 



, W i<n, (74) 

,3i = l,jji,=0, Vfe<i,t;fc = l,Vfe>i 



where we have denoted the objective function by J-psmto- It is clear form these values that, for (i — I) < rg < i, 
the steepest descent is along decreasing Vi with /3i = 1 and 71 = a. Putting this in equation ( |72] | we get 

n 

l'^^) =min5^(4n + 1 - 2i)v, + n\a - 1); (75a) 

1=1 

n 

constrained by: — Vi)^ < rs\ (75b) 

i=l 

< Ul < • • • < Un. (75c) 

Now, using Lemma [3] we obtain the minimum value of the the above optimization problem, which can be written 

as 

dl%°n,n^n){rs) = rfn,3n(r.) + (a - 1 ) , < r. < n. (76) 

Next, we obtain the optimal value of the objective function for values of rg > n. Note that for any 7'^ > n, 
Vi = {) y i, which along with equation ( I72fb imply that (5j = 1, -fj > 1 Vj. Putting this in equation ( |72] | we get 

n 

dl%:n,n,n) ('^.) = mill ^(2n + 1 - 2^)J, - n\ (77a) 

i=l 

n 

constrained by: — < {rg — n); (77b) 

i=l 

1 < 71 < • ■ ■ < 7n. (77c) 

To bring the above problem into a form amenable to Lemma [3] we use the following variable transform in the above 
set of equations: 7^- = 7,: — 1. This results in the following equivalent optimization problem. 

n 

dl%:n,u,n)irs) = min^(2n + 1 - 2z)7,', (78a) 

n 

constrained by: — 1 — 7-)^ < {r^ — n); (78b) 

i=l 

0<7; <---<7;, (78c) 
which in turn by Lemma [3] attains the following optimal value: 

= (a - l)d„^„ ( ITI^Z^ ) , n<rs<na. (80) 



{a -I) 

Finally, combining equations ( |76] l and dSTj i we obtain the desired result 
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Appendix D 
Proof of Theorem|4] 



Substituting Mi = AI2 = M and = 1 in equation (|26] | we obtain 

M M M 

df(M,N,M.N,)(rs) = ^(2M + TVi + iV2 + 1 - 2t)u, + ^(M + iVi + 1 - 2j)/3, + ^(M + TVa + 1 - 2fc)7fc 

i=l j=l k=l 

M (M-k) M min{(Ni-j),M} 

~{M + N2)M + Y, E (l-«.-7fc)++E E {l-v^~|3,)+■, (81a) 

fc=l 4=1 i = l 4=1 

P 91 92 

constrained by: E("2i - v.,)'^ + E("ii ~ ft )^ + E('^22 - Jk)'^ < (81b) 

i=l j=l fc=l 

0<ui<---<Up; (81c) 

0</3i <••• (81d) 

0<7i < ••• <792; (81e) 
+ /?,)> 1, V(z+j) > (iVi + 1); (81f) 

(^^4 + 7fe) > 1, V(i + fc) > (Af + 1). (81g) 

Differentiating the objective function with respect to Uj, V i,/3i and 71 we find that 



a-FCSIT 



9Ui 



/3i = l,7i = l 



-FCSIT 



7i=l,afc=0, Vk<i,ak = iyk>i \f i < n, (82) 
;9i=l,Qfc=0, VA:<i,ak=l,Vfe>i ' 



where we have denoted the objective function by J^a-FCSiT- It is clear form these values that, for (i — 1) < < z, 
the steepest descent is along decreasing Vi with /3i = 71 = 1 and i < n. Putting this in equation (ISTT i we get 

n 

C(A/Vi,M,iv.) i^s) = min E(2A^ + N1+N2 + I- 2i)v, + M(7Vi ~ Af); (83a) 

1=1 

j\/ 

constrained by; ^ ^«)^ — ''s' (83b) 

i=i 

< wi < • • • < w„. (83c) 

Now, using Lemma |3] we obtain the minimum value of the the above optimization problem, which can be written 

as 

<X*^iVi.M,JV.)('^s) = duM+N^+N^rs) + M{Ni - A/), < r, < M. (84) 
Next, we evaluate the optimal value of the objective function for values of rg > AI. Note that for any Vg > Af, 
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= V i, which along with equation (81gi imply that 7^ > 1 Vj. Putting this in equation dSH i we get 



dftAlN,,Ai,N,)(.^^) =min^(A/ + TVi + 1 - 2j)/3, - (85a) 



constrained by: ^^(a — 7i)^ ^ (''s — "■); (85b) 



4=1 



1 < 71 < • • • < 7n. (85c) 

To bring the above problem into a form amenable to Lemma [3] we use the following variable transform in the above 
set of equations: 7^' ~ ji — 1. This results in the following equivalent optimization problem. 



<!Kn,„,„) (rs) = min ^(2n + 1 - 2iH, (86a) 
1=1 

n 

constrained by: — 1 — 7j')^ < ('^s ~ J^); (86b) 

!=1 

< 7; < . . . < 7;, (86c) 
which in turn by Lemma [3] attains the following optimal value: 

dl%:n,n,n)irs) " (^3^) - < < n (87) 

= (a ~ l)d„,„ -r^ — — ] , n < Ts < na. (88) 



(a-1) 

Finally, combining equations (|76] | and dSTl i we obtain the desired result. 

Appendix E 
Proof of Lemma |5] 

Let us denote the event that a target rate tuple (ri log(p), r2 log(p)) does not belong to 7?.iml by Oiml, i-e., 

{•H : (ri log(p), r2 log(p)) i 7^IML}. 

Now, let us denote the maximum among the average probability of errors at both the receivers be denoted by Ve, 
then using Bayes' rule we get 

Ve =^e|o,ML MOiml} + V.^o^ ^^{O^^m} , (89) 

<Pr{OiML}+^e|o,w, (90) 

where Ve\£ denote the conditional average probability of error given the event £. Note that, the above equation 
holds for any SNR. When the target rate tuple belongs to TZiml, letting the block length be sufficiently large the 
probability of error given Of^L be be made arbitrarily close to zero. Therefore, letting the block length of the 
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code goes to infinity at both side of the above equation we obtain 

< Pr{OiML} =Pr {{Ic, < n log(p)} U < r2 log(p)} U {I,, < log(p)}} , (91) 
< Pr{{/,. <raog(p)}}, (92) 

i=l,2,s 

(a) 

= max Pr < n log(p)}} , (93) 

i— l,2,s 

= max p~''™'i-fl,iVi,A/2.«2)(''') = p~™'"'=1.2,s{4"(Il/i,]Vi,Af2.iV2)(''i)}^ (94) 

i— l,2,s 

where step (a) follows from the fact that in the asymptotic SNR the largest term dominates and the last step follows 
from equation ([32l). Finally, the desired result follows from the fact that 7'e=p^'*(M.«i>J>^-"2)('"i'''2). 

Appendix F 
Proof of Lemma[6] 

The joint distribution of {^,v) can be obtained from equation dSOb substituting fw^lWsi-) = 1- The rest of the 
proof follows the same steps as the proof of Lemma |4] and hence skipped to avoid repeating. 

Appendix G 
Proof of Theorem[6] 

Since the FCSIT DMT is an upper bound to the No-CSIT DMT, it is sufficient to prove that when the condition 
of equation dTIT i is satisfied the expressions of Theorem |2] and Lemma [T] are identical. 

It is clear from the comparison of equations ( |29] l and ( l40l i that for = 2r > n they are identical for all values of 
a > 1 and hence when a satisfies equation (|4TJ. Now, it is only necessary to find a condition when the expressions 
in equations ( |29] l and ( |40] | are identical even when < n, which is what we derive next and turns out to be 
identical to the condition of equation (|4TI ). 

It is clear from equation ( |40b and ( |29] ) that, 

rf™n,n,„,n)(^.) < dlf^l.,n,n)irs), for r, < n. 

Therefore, the DMTs given by Theorem|2]and Lemma|2]are identical only if for < n, the single user performance 
is dominating, i.e., 

dn,n{r) <c?™(n,„^„,,i)(^s = 2r); 

duAr) <d«,2n(2r) +n^{a~l), (95) 

where in the last step we have substituted the value of rf™^ „ „ from equation ( l40l ). Since c?™^ „ „ n)(2'') 

decays much faster than cJ„_„(r) with increasing r and both are continuous functions of r, equation ( |95] l will be 
valid for all r < 2. if it is true for = §. Substituting this in equation ( |95] l we get 

dn,n <n^(a - 1), 
or, a >1 H TT-^. 
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